Tunable optical Aharonov-Bohm effect in a semiconductor quantum ring 
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By applying an electric field perpendicular to a semiconductor quantum ring we show that it is 
possible to modify the single particle wave function between quantum dot (QD)-like to ring-like. 
The constraints on the geometrical parameters of the quantum ring to realize such a transition are 
derived. With such a perpendicular electric field we are able to tune the Aharanov-Bohm (AB) effect 
for both single particles and for excitons. The tunability is in both the strength of the AB-effect as 
well as in its periodicity. We also investigate the strain induce potential inside the self assembled 
quantum ring and the effect of the strain on the AB effect. 
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I. INTRODUCTION 

The electron and hole wave functions acquire an extra 
phase when moving in the presence of a perpendicular 
magnetic flux. The acquired extra phase will be different 
for electrons and holes, and this phase difference can be 
observed through photoluminescence (PL) experiments, 
which exhibits an optial Aharonov-Bohm (AB) effect. 
Experimental verification of this effect has been reported 
in PL measurements of radially polarized neutral excitons 
in a type-II quantum dot structure^!. However, the op- 
tical AB effects will be strongly suppressecP when both 
the electron and the hole are spatially confined within 
the same region due to the Coulomb interaction, i.e. in 
a quantum ring (type-I quantum dot). This has spurred 
a considerable search and study of the optical AB effect 
of neutral excitons in semiconducto r qu antum rings by 
both theoreticaP^ and experimental^^ groups. 

The optical AB effect in a quantum ring can be en- 
hanced when the exciton is radially polarized, either by 
the application of an external electric fielcP^ or due to 
a radial asymmetry in the effective confi neme nt for elec- 
trons and holes^" 8 -'. Theoretical studie d 3 ^ 8 ! 13 ! on one di- 
mensional rings predict that for a quantum ring whose ra- 
dial size is comparable to the exciton Bohr radius (i.e. in 
the weakly bound regime), the ground state energy could 
display nonvanishing AB oscillations. However, numeri- 
cal calculations on two dimension narrow rings^ showed 
that: 1) there is no observable AB oscillation in the ex- 
citon ground state energy, but 2) AB oscillation can be 
present in some low-lying excited energy levels!^. 

Here we study a nano-ring structure, with and without 
strain, which has a varying height in the radial di rection 
as is often found in self-assembled quantum ring a 12 ! 1 I 15 l 
Applying a strong electric field in the perpendicular di- 
rection polarizes the exciton. We will show that the AB 
effect is strongly affected by the shape of the nano-ring 
structure and the external electric field which has a large 
effect on the Coulomb interaction between the electron 
and the hole, especially for quantum rings with a large 
height. This gives us the possibility to tune not only the 
shape of the electron and hole wave function, but even 




FIG. 1: (Color online) Schematics of the volcano-like quan- 
tum ring. A perpendicular electric and magnetic field are 
applied in the vertical direction. Bottom: the structure of 
the quantum ring in the p-z plane. 



the magnetic field induced oscillations in the exciton en- 
ergy through the application of an external electric field. 

This paper is organized as follows: we present the phys- 
ical model in Sec. [TTJ In the first part of Sec. |III| we use 
the finite element method to calculate the electron and 
the hole wave functions and energies in unstrained semi- 
conductor quantum rings. Subsequently we calculate the 
total energy of the exciton in the presence of Coulomb 
interaction by diagonalizing the total Hamiltonian in the 
space spanned by the product of the single particle states. 
We show the tunability of the AB effect for both the sin- 
gle particles and the exciton energy by applying a per- 
pendicular external electric field. In Sec. |IV[ we will give 
the equivalent results for a strained semiconductor ring, 
and we will show how the strain affects the results. Our 
conclusions are presented in Sec. V. 



II. MODEL 

Figure. [T] shows the investigated geometry of the three 
dimensional ring, A volcano shaped ring is embedded in a 
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barrier material which is different from the ring material. 
In a three dimensional semiconductor quantum ring, the 
full Hamiltonian of the exciton within the effective mass 
approximation is given by 

Htot = { P i - lity 2^7 ( p i - + Vc - fh) 

j=e,h 3 

+ ]T SE j(rj) + E v i<fi) - eFz e + eFz h , (1) 



j=e,h 



j=e,h 



where Vj(f f j) is the confinement potential of the elec- 
tron (hole) due to the band offset of the two materi- 
als which will be different inside and outside the ring. 
V c (r e — r h ) = e 2 /4ire\ (f e — r h )\ is the Coulomb poten- 
tial between the electron and the hole, and SEj(fj) is the 
strain-induced shift of the energy which depends on the 
strain tensor e^ . We did not take the piezoelectric po- 
tential into account since it is negligible compared to the 
other terms in our case. The last two terms of Eq. (!]) are 
the potential energy in the presence of the perpendicular 
top to bottom directed electric field F. 



only the lowest electronic subband and the highest hole 
band (heavy hole) is occupied. In GaAs, the electron 
and the hole have effective mass m e /mo = 0.063 and 
trik/mo = 0.51, respectively. The static dielectric con- 
stant is £ = 12.5eo and the band gap is E g = 1.42 
eV at helium temperature. While for AlGaAs, we have 
nT-e/ m o — 0.082, m^/m = 0.568, e = 12.5eo, an d a band 
gap of Eg — 1.78 eV. This results in a band gap difference 
of AE g — 360 meV between GaAs and AlGaAs, which 
leads to a conduction band offset of about AE C = 250 
meV and a valence band offset of about AE V = 110 meV, 
so we can take V{f e (h\) — inside the ring and AE C ^ 
outside the ring. The parameters used were taken from 
RefP^. As the dielectric constant is practically the same 
everywhere and the difference of the lattice constant for 
GaAs and AlGaAs is very small, we may ignore the di- 
electric mismatch effect, and the strain induced term in 
Eq. Q. 



III. GaAs /AlGaAs QUANTUM RING 

A. Single particle energy and wavefunction 

As a model system we consider first a volcano shaped 
GaAs ring surrounded by AlGaAs. We assume that 



First, we calculate the single particle energy and the 
corresponding wave function. Because of cylindrical sym- 
metry we rewrite the Hamiltonian in cylindrical coordi- 
nates. Assume the wave function of the single particles 
to be #(p, z, 9) — Tp e (h){p, z)e~ ll < h '> 8 , after averaging out 
the angular part of the wave function, we obtain the 2D 
single particle Schrodinger equation : 
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+ V(f e{h) ) +qFz e(h) )ip e ( h )(p,z) = E PtZ ip eW (p,z), (2) 



here q is — e for the electron and e for the hole. Tak- 
ing advantage of cylindrical symmetry we only need to 
solve the problem in a half section of the ring as shown 
at the bottom inset of Fig. [l] We use the finite element 
methods to obtain the ground state energy for different 



angular moment l e (h) an d electric field F as a function 
of magnetic field B. For a better understanding and fur- 
ther consideration of the transition, we first write down 
the single particle Schrodinger equation in dimensionless 
form (assume B = and l e (h)=o)'- 
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) + V'(r e{h) ) +A F F'z e (t l) j i> e (h)(p,z) = E' pz ijj e{h) (p,z), 



(3) 



where E', = E, 



R 2 is the 



- P ,z - ~p,z/E , and E = h 2 /2m eo(ho 
energy unit. We take Rq = 1 nm as the length unit, 
and m eo (/j ) is the effective mass of the electron (hole) in 
GaAs. The effective band offset V'(r e ( h )) = V(r e ( h ))/E 
is 0.425 for the electron and 1.472 for the heavy hole. We 



take the unit of the electric field to be Fi = 1 kV/cm and 
the coefficient A F — FiR e/E a = 0.00263m eQ (; lo ). 

Fig. [2] shows the single particle energies of the electron 
and the hole as a function of the perpendicular magnetic 
field B for different values of the angular momentum and 
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FIG. 2: (Color online) Electron and hole energy for different 
electric field F and angular moment l e (h) as a function of 
magnetic field. The size of the quantum ring is hi = 4 nm, 
tl2 = 6 nm, _Ri = 12 nm and R2 = 30 nm. The arrows 
indicate angular momentum transitions in the single particle 
ground state. 



the electric field F. Here the size of the ring is chosen to 
be hi = 4 nm, h 2 = 6 nm, R\ = 12 nm and R 2 = 30 nm. 
Notice that the absolute value of the angular momentum 
of the ground state for both the electron and the hole in- 
creases with B and a pronounced AB oscillation is found 
(the second transition of the hole does not take place be- 
low B = 30T, but it will definitely be there above 30T). 
It is clearly seen that the electric field has a larger effect 
on the hole than on the electron. Furthermore, the elec- 
tric field pushes the electron and the hole wavefunction 
in opposite z-direction which results in a difference in the 
shape of the wave functions. When we decrease the elec- 
tric field or even change the direction of the electric field, 
the angular momentum transition of the electron ground 
state will be more obvious, while for the hole, the transi- 
tion becomes invisible (when F = 100 kV/cm, there is no 
transition for the hole ground state for a magnetic field B 
below 30 T). This will provide us with the possibility to 
tune the Aharonov-Bohm effect using the perpendicular 
electric field F, i.e. we can turn the AB effect off and 
on with this perpendicular electric field. We can have 
a ring-like wave function (i.e. the probability is zero in 
the center of the ring) or a (QD-like wave function (i.e. 
not zero probability of the wave function in the center) 
depending on the perpendicular electric field F. 

Fig. [3] shows the contour plot of the electron and the 
hole wave function in the p — z plane. When the electric 
field changes from 100 kV/cm to —100 kV/cm the wave 





FIG. 3: (Color online) Contour plot of the ground state wave 
function of the electron and the hole for different values of 
the perpendicular electric field F in the (p, z) plane. Here 
B — 0, l e th) = and the electric field F from top to bottom 
are 100 kV/cm, 50 kV/cm, 0, -50 kV/cm, -100 kV/cm. The 
contour of the ring is shown by the black curve. 



function of the hole changes from a QD-like to a ring-like 
wave function, while the wave function of the electron 
changes from ring-like to QD-\ike. Notice that the wave 
function of the hole is more sensitive to the electric field, 
which can be understood from the dimensionless form of 
the Schrodinger equation, Eq. Notice that the term 
related to the electric field is 0.00263m eo (/ lo )F'z, which 
is proportional to the effective mass. As the heavy hole 
has a larger effective mass, the effect of the electric field 
will be larger. 

From a careful inspection of Eq. |3| and Fig. [3] we 
determine the constraints in order to have an electric 
field tunable Aharonov-Bohm effect: 1) The electric 
field term ApF'dz must be comparable to the differ- 
ence of the confinement energy of the QD-like and ring- 
like wave function, which is proportional to l/h\ — l/h%. 
When passing from the Q-D-like to the ring-like wave 
function, the difference of the z coordinate is propor- 
tional to {hi — hi)/2, so the first constraint condition is 
ApF'(h 2 — h\)/2 ~ 1/hf — \jh\- For example in our 
case, hi = 4 nm, h 2 — 6 nm and for the heavy hole 
with an effective mass = 0.51too the electric field 
should be at least 2(l//if - l/hl)/A F /{h 2 - hi) * F t = 
2(hi+h 2 ) /(hih^Ap)* Fi, which is several tens of kV/cm. 
If we decrease hi, then the electric field should be larger. 
We see from Fig. Qa) that if hi = 2 nm, the wave func- 
tion of the hole is still ring-like when F — 100 kV/cm, 
and we need to increase the electric field in order to have 
a QD-like wave function, which leads to the second con- 
straint condition: 2) the confinement energy of the elec- 
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FIG. 4: (Color online) Contour plot of the ground state wave 
function of the heavy hole for different sizes of the ring. The 
size (in unit of nm) and the corresponding electric field are 
specified above each figure. 

tron (hole) when it is QD-like (i.e. which is ~ V 1 ?) 
should be smaller than the band offset V'(f e (/ l )) of the 
quantum ring, otherwise the electron (hole) will move out 
of the ring to the barrier for large electric field values, as 
can been seen from Figs. |4jb,c). The wave function is 
still ring like when F — 200 kV/cm, but when we in- 
crease the electric field to F = 250 kV/cm, the wave 
function starts to penetrate in the barrier region. The 
AB effect is still tunable, but the effect of the electric 
field will be very small in case of a small ring height hi- 
3) The confinement energy difference of the QD-like and 
ring — like states, which are proportional to \jh\ — \jh\ 
should not be very small, otherwise the wave function of 
the electron (hole) will always be QD-like and shows no 
observable AB effect (as we can see from Fig. |4]jd), when 
hi = 5.5 nm and h 2 = 6 nm the wave function is still 
QD-like when F = —100 kV/cm). 4) In order to have a 
QD-like wave function, the radius of the ring R\ should 
not be very large compared to the height of the ring (oth- 
erwise the electric field will not have a significant effect 
on the wave function), and this constraint is much more 
important when the height of the ring is small (we can 
see this by comparing Figs. [IJe, f) with Figs. |4ja) and 
[3]). The above four constraints were obtained for holes, 
but they are also valid for electrons. 



B. Exciton energy including Coulomb interaction 

The total exciton energy is calculated using the con- 
figuration interaction (CI) method. We first construct 
the total exciton wave function as a product of linear 
combinations of single electron and hole wave functions, 
and then calculate the matrix elements of the exciton 
Hamiltonian. By diagonalizing the obtained matrix we 
obtain the exciton energies and their corresponding wave 
functions. 

The Hamiltonian of the exciton is 

H tot =H e + H h + U c , (4) 
here H e and Hh are the single particle Hamiltonians, 

U c = U * . 1 =, 

is the Coulomb energy between the electron and the hole, 
where Uq — — e 2 /AiteRq. 

As the total angular momentum is a good quantum 
number because of cylindrical symmetry, we assume the 
wave function of the exciton to be: 

*L (r e , f h ) = J2 C ^k (f e) r h ) , (6) 
k 

for given total angular moment L. Here, (r e , r),,) = 
i/jn e e~ d '' (p '' ipn h e~ llhiph , and the exciton wavefunction is 
constructed out of single particle eigenstates, and k 
stands for the indices (n e ,n/j,Z e ,//s). l e + lh = L should 
always be satisfied for fixed L, n e (rih) is the quantum 
number of the single particle radial wave function. As 
the energy of the single particle eigenstates with large 
angular moment quantum number l e (//,) and quantum 
number n e (nh) is much larger as compared to the ones 
with smaller angular moment quantum number, we can 
limit ourselves to several tens of single particle eigen- 
states. With this wave function, we can construct the 
matrix of the total Hamiltonian and by diagonalizing the 
obtained matrix, we can get the eigenvalues and eigen- 
vectors. 

For a fixed total angular moment, we have 



U c (k,j) = C k C 3 (<f> k (r e ,r h )\U c \<f> j (f e ,r h )) = C k C 3 J J J J ip n ^n h ipm e ^m h Ai eJ j h jj ek ,i hh p' e dp' e p h dp' h dz' e dz' h , (7) 
and where Em is the angular part of the integral 

r 2 ^ r 2 ^ e -i(lej+lhj-lek-lhk)'Pe e i(lhj-lhk)(<Pe-<Ph) 

Ai ej ,i hj ,i Bk ,i hk = Uq I I dtp e d<p h —j= = -- =. (8) 

V + Ph - 2 P'eP'h COS (^ - <Ph) + (K - A) 
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FIG. 5: Exciton energy for total angular momentum L = 
as obtained from the CI method (solid line) and compared 
to the results of the ground exciton energy without Coulomb 
interaction (dot-dashed line), and the ground state energy by 
taking the Coulomb energy as a perturbation (dashed line). 
Here the top figure is for F — 0, and the bottom one is for 
F = 150 kV/cm. 



We know from Eq. ([8]) that, in order to make the inte- 
gral nonzero, l e j + Ihj — l e k ~ hik should be zero, which 
indicates that the total angular moment is a good quan- 
tum number. If we take tp e — iph as a new variable in 
Eq. Q , it can be changed into an elliptic integral, which 
is easy to calculate numerically. The remaining inte- 
gral is easy to calculate numerically. The effective ex- 
citon Bohr radius is = Airsh 2 / (fie 2 ) — 11.7673 nm 
(fi = l/(l/m e + 1/m/j)), and we take a small volcano 
shaped ring (r c = 10 nm, R = 16 nm, h\ = 2 nm and 
hi = 4 nm), whose size is comparable to the exciton 
Bohr radius. The results of the total exciton energy for 
F = and F — 150 kV/cm are shown in Fig. [5j and 
as a comparison, we also show the results of the ground 
exciton energy without Coulomb interaction and the to- 
tal ground state energy by taking the Coulomb energy 
as a perturbation (we took only the lowest electron and 
hole state into account, the result is much closer to the 
real one for a smaller ring when the Coulomb energy is 
small as compared to the kinetic energy) . From Fig. [5] 
we know that the Coulomb interaction has a large effect 
on the total exciton energy, and unlike the perturbation 
result, the exciton energy from the CI method shows a 
monotonous increase with magnetic field. There are no 
oscillations as function of the magnetic field. The exis- 
tence of the top to bottom directed electric field decreases 
the exciton energy by tens of meV but does not bring any 
oscillation. 
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FIG. 6: (Color online) Ground state energy and the five lowest 
excited exciton energies as a function of the magnetic field. 
The top figure is for the case of F = 0, and the middle one 
is for F — 150 kV/cm. (c) Second derivative of the exciton 
energy with respect to the magnetic field. 



Figures [6]Ja, b) show the CI results, for the exciton 
energies with different values of the total angular mo- 
mentum L, when F = and F = 150 keV. It shows that 
the state with total angular moment L = is always the 
ground state, which is different from the perturbation 
theory result shown in Fig. [5] where the total angular 
moment of the ground state switches between L = and 
L = 1. We notice that the energy of the exciton states 
show no oscillation with increasing magnetic field, both 
with or without perpendicular electric field. The effect of 
the electric field on the exciton leads to a decrease of the 
exciton energy and a larger difference between the energy 
of the lowest exciton states. In Fig. |6jc) we plot the sec- 
ond derivative of the exciton energy with respect to the 
magnetic field, for both with (solid curves) and without 
(dotted curves) electric field. We find that the second 
derivative of the exciton energy shows a clear oscillation 
with respect to the magnetic field, both for the ground 
state and excited states energy. But since the size of the 
ring, as we specified above, is very small, the period of 
the oscillation is large. The oscillation in d 2 E to t/dB 2 is 
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FIG. 7: (Color online) Exciton ground state energy as a func- 
tion of the magnetic field for F = 0, F = —200 kV/cm and 
F — 200 kV/cm (inset). The curves in different color stand 
for different y scale. 



slightly more pronounced when a perpendicular electric 
field is applied. The reason is that the electric field makes 
the electron and the hole move in the opposite directions 
which enhances the polarity of the exciton. 

Till now, we did not observe any optical AB effect in 
the exciton ground state energy and any large influence 
of the perpendicular electric field on the oscillation. In 
order to see the effect of the electric field more clearly, 
we calculate the energy for a quantum ring with a larger 
height and a much larger electric field. Since the pe- 
riod of the oscillation for the previous case is large, we 
increase the radius of the ring in order to decrease the 
period. Fig. [7]Ja) shows the exciton ground state energy 
with respect to the magnetic field in the case of F = 
kV/cm, F = —200 kV/cm and F = 200 kV/cm (inset of 
Fig. [7]). Here R\ — 14 nm, R2 = 18 nm, hi = 1 nm and 
h 2 = 10 nm (the size is still comparable to the effective 
Bohr radius). As the state with total angular moment 
L = always has the lowest energy, the ground state is 
the same as the L — state here . It is clear that al- 
though the F — state shows no oscillation at all, the 
ground state energy in the presence of an electric field 
shows a weak but observable oscillation, especially for 
the case of F = —200 kV/cm. And the ground state en- 
ergy is smaller in the presence of a top to bottom directed 
electric field. 

For a better understanding of the existence of the 
ground state AB oscillations and the difference of the 
oscillation period for different values of the electric field, 
we calculate the average value < p > in the radial di- 
rection and < z > in the z direction of the electron and 
the hole in the exciton ground state, which we show in 
Fig. H From Fig. ^a) we see that < p > and < z > for 
both electron and hole have no clear step-like behavior, 
but they rather change smoothly. Furthermore, the dif- 
ference of the average values of the electron and the hole 
in both radial and z direction is rather small, especially 
in the z direction. This makes then the Coulomb inter- 
action energy very large and the polarity of the exciton 
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FIG. 8: (Color online) Average distance < p > (left y scale 
in same color) and < z > (right y scale) in unit of nm of the 
electron (hole) as a function of magnetic field for (a) F = 0, 

(b) F — —200 kV/cm and (c) F = 200 kV/cm. In (b) and 

(c) the curves of the electron and the hole are separated in 
different sub figures. 



extremely small, which greatly weakens the oscillations 
in the exciton ground state energy. That is the reason 
why we can not observe any oscillation of the exciton 
energy in the absence of the electric field. But in the 
presence of a perpendicular electric field, we see from 
Figs. [8^b) and [8|c) that the difference of the average 
values between the electron and the hole becomes much 
larger. And the average values < p > and < z > show 
pronounced increasing (or decreasing) step-like behavior 
with increasing magnetic field, especially for the electron. 
In both cases of F = 200 kV/cm and F = —200 kV/cm, 
< p e > and < z e > show obviously oscillations with 
increasing magnetic field, and their oscillation period is 
more or less the same as the period of the oscillation 
in the exciton energy. That is because the coupling of 
the electron and hole states with different angular mo- 
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mentum becomes much smaller as a result of the weak 
coulomb interaction. If we look at the normalization con- 
stant Cfc in Eq. ^ we know that the step-like behavior 
in < p e > and < z e > in Figs. [8jb,c) originate from the 
angular momentum transition of the main contributing 
basis function in the total exciton wave function. With 
increasing magnetic field, the angular momentum pair 
(l e ,lh) in the state $fc (r e ,fh) which has the largest con- 
tribution to the total wave function changes from (0, 0) 
to (—1,1) or to even larger values, and this transition 
becomes more notable in the presence of a perpendic- 
ular electric field as the coupling between the different 
contributing basis functions is much weaker. 

From Fig. [8] we also see that there are some differ- 
ences between the cases F — —200 kV/cm and F = 200 
kV/cm. In the case of F = —200 kV/cm, the electron 
is in the center area of the ring while the wave func- 
tion of the hole has a ring shape. The Coulomb inter- 
action makes the electron (hole) move to the top (cen- 
ter) area of the ring towards the hole (electron) . When 
the magnetic field is increased, the electron (hole) wave 
functions change towards the area with larger (smaller) 
values of p in which the hole (electron) is located, but 
as a result of the strong electric field and confinement 
energy, only < p e > is mainly changed. Moreover, com- 
pared to the electron, the averages < p > and < z > 
for the hole change more smoothly with increasing mag- 
netic field. The behavior for F = 200 kV/cm is different. 
The electron now is in the top area of the ring, p e and 
z e change over a rather small range as a result of the 
strong confinement in the top area of the ring. More- 
over, with increasing magnetic field, < p > and < z > 
for both the electron and the hole increase towards the 
area with larger value of p, but the change for the electron 
is very small. This is different from the case of F = —200 
kV/cm. 

In Fig. |] we show the second derivative of the exciton 
ground state energy with respect to the magnetic field 
for different values of the perpendicular field. We clearly 
notice that the oscillation is enhanced in the presence of 
the perpendicular electric field, independent of the direc- 
tion of the electric field. The oscillation is much larger 
when the electric field is directed from bottom to top. 
This can be understood by comparing Figs.[8][b,c), where 
the < p > (< z >) differences between the electron and 
the hole are both larger. But the difference is smaller 
in the case of F = 200 kV/cm, where the Coulomb en- 
ergy U c is larger, which weakens the AB oscillation in 
the exciton ground state energy as compared to the case 
of F = -200 kV/cm. We also verified in Fig. [9] that 
the period of the oscillation is larger when increasing the 
bottom to top directed electric field, but smaller by in- 
creasing the top to bottom directed electric field. The 
period change is a consequence of the change of average 
radius of the electron and/or hole with changing electric 
field. However, it is hard to obtain this result only from 
the changing in < p e > and < p) % > alone. As we notice 
from Figs. J8^b) and [8^c) that < p e > and < z e > is 
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FIG. 9: (Color online) Second derivative of the exciton ground 
state energy with respect to the magnetic field for different 
values of the perpendicular field F — (black solid line), 
F = -100 kV/cm (red solid line), F = 100 kV/cm (red dash- 
dotted line), F = -200 kV/cm (light blue solid line), F = 200 
kV/cm (light blue dashed line) and F — 250 kV/cm (dark 
blue dotted line) 



smaller in the case of F = —200 kV/cm as compared to 
F = kV/cm, but < ph > and < Zh > is larger, which 
probably decreases the AB oscillation period of the total 
exciton energy. Moreover, although < p e > and < z e > 
is larger in the case of F = 200 kV/cm, but < ph > 
and < Zh > is smaller. The point is that as the hole 
has a much larger effective mass, the average radius of 
the exciton is mainly determined by the electron. As we 
explained in Appendix. B, the period of the oscillation 
is proportional to m%l < p e > 2 +m e / < ph > 2 , and is 
mainly determined by the effective radius of the electron. 

In order to study the exciton recombination, we cal- 
culate the oscillator strength (the dimensionless quantity 
that expresses the strength of the transition) for the state 
with total angular momentum L = 0. The os cillato r 
strength for the exciton ground state is defined a d 16 * 17 * 



fy 



2 | < ex|]T>„|0 > 



E, 



Eq 



(9) 



where m is the free electron mass and \ex > (|0 >) is the 
exciton state (single electron and hole pair). By using 
the envelope-function approximation, we can derive^ 



2P 2 

m(E ex - E ) 



* s {re,r e )dr c 



(10) 



here P includes all intra matrix-element effects, 
(^h, r e ) is the exciton ground state wave function. For 
simplicity, here we will focus on the main variable part of 
the oscillator strength which is Of = \J 4" g (f e ,r e )dr e \ , 
named the overlap integral. The result of the overlap in- 
tegral for the exciton ground state is shown in Fig. [10} for 
F = -200 kV/cm, F = kV/cm and F = 200 kV/cm. 
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FIG. 10: (Color online) Overlap integral Of = 
|J V& 9 (r e , r e ) dr e | 2 of the exciton ground state energy as a 
function of the perpendicular magnetic field. The pink dashed 
line is for F = 0, red solid line stand for F = —200 kV/cm 
and black dash-dotted line is for F = 200 kV/cm. Note the 
different scales for the three curves 



In the presence of the electric field, Of exhibits more 
structure, and we can clearly see steps in the overlap in- 
tegral. The steps in the overlap integral have the same 
period as the AB oscillation, which also originates from 
the angular momentum transition. Moreover, the over- 
lap integral for the case of F = is much larger, which is 
a consequence of the extremely small polarity of the ex- 
citon. Of in case of F = —200 kV/cm shows the smallest 
value as compared to the other two cases, which means 
the exciton has the largest polarity and is most stable, 
and the AB effect is the strongest (which confirm our re- 
sults of Fig. [9]). We also find that for all the three cases, 
the overlap integral increases with increasing magnetic 
field, which means the exciton has a smaller polarity and 
also a weaker AB oscillation. This is quite different from 
the case of a single particle in a one dimensional ring 
where the AB oscillation has almost the same strength. 
The reason is that by increasing the magnetic field, the 
Zeeman term makes the angular momentum (absolute 
value) of both the electron and the hole increase, and 
consequently the electron and the hole will have a larger 
effective radius (note that in the one dimensional ring 
it is always fixed). However, because the confinement 
potential is different inside the 3D ring it prevents the 
electron and the hole to move to the outer part of the 
ring which has a radius larger than the top area of the 
ring. Which leads, as shown in Fig. [8j to a smaller dif- 
ference in the electron and hole position inside the ring 
with increasing magnetic field, and consequently, a larger 
value of the overlap integral and weaker AB oscillation. 
We can also verify this result in Fig. |9j where we see that 
the amplitude of the AB oscillation decreases when we 
increase the magnetic field. 



IV. InGaAS/GaAs QUANTUM RING 

Now let us turn to a different system where strain is 
important: Ini_a;Ga x As ring surrounded by GaAs, the 
concentration of Ga is proportional to the coordinate z 
inside the ring!--, which is x = 0.4 — 0.05z. We take 
R\ = 15 nm, R2 = 22 nm, hi = 0.5 nm and I12 = 4 nm 
for the ring. 

For Ini^GazAsP, we have the effective masses 
m e /m = 0.023 + 0.037x+0.003a: 2 , m h /m = 0.41+O.lz, 
dielectric constant e = (15.1 — 2.87a; + 0.67x 2 )eo, and a 
band gap of E g = 0.36 + 0.63a; + 0.43a; 2 eV. This results 
in a band gap difference of A£" g = 1.06 — 0.63a: — 0.43a; 2 
eV between GaAs and Ini-^Ga^As, we take 25% of AE g 
be the valence band offset and 75% be the conduction 
band offset. Since x = 0.4 — 0.05z inside the ring, the 
band gap difference will be the largest at the top of the 
ring (which is (1.06 - 0.63 * (0.2) - 0.43(0.2) 2 ) = 0.9168 
eV), we assume V(fi l ) = when z = 4 inside the ring, 
then we find that V{f h ) = 0.25 * 0.9168 = 0.275 eV out- 
side the ring and V(r h ) = 0.275 - 0.25(1.06 - 0.63a; - 
0.43a; 2 ) = 0.05324 - 0.01461z + 0.00032z 2 eV inside 
the ring, while for the conduction band offset we have 
V(f e ) = 0.5032 + 0.9168 * 0.75 = 1.145 eV outside the 
ring and V(r e ) = 1.145 - 0.75(1.06 - 0.63a; - 0.43a; 2 ) = 
0.62756 - 0.03409z + 0.00075z 2 eV inside the ring. For 
convenient of our calculation, we take the value of the 
hole confinement potential V(fh), as in case of the elec- 
tron, which make the hole also be confined in a potential 
well, larger value of V(fh) stands for higher energy. As a 
simplification, we do not take the dielectric mismatch ef- 
fect into account because it is nearly the same in the ring 
and in the barrier material, but just assume e — 12.5eo 
inside and outside the quantum ring structure. 

As in Sec. |III| we will solve the single particle 
Schrodinger equation first. As the lattice constant!^ is 
ai = 0.56533 nm for GaAs and a 2 = 0.60583 - 0.0405a; 
nm for In^^Ga^As (x from 0.2 to 0.4 ), there is a lat- 
tice mismatch (02 — a\)/ai of 6%, which results in a 
large strain. The difference from the previous case is 
that in a strained quantum ring the total confinement 
potential terms V(r e ^) in Eq. (|2j) now comes from the 
band offset energy due to the band gap difference, and 
the strain induced term. We calculated the strain by 
adapting a method developed by John Davi es which is 
based on Eshelby's theory of inclusion d 19 * 20 *, where the 
elastic properties are assumed to be isotropic and ho- 
mogeneous (for more details, see Appendix. C). In our 
model the lattice mismatch between the two materials 
is £ = (0.243 + 0.02025z)/5.6533 inside the ring and 
outside. By using the finite element method, it is easy to 
obtain the eigenstrain. The strain in our case will change 
the potential of the electron (hole) thus modify the band 
structure including the conduction and the valence edge 
energies which are among the most important parameters 
characterizing them. 

We assume that the conduction band is decoupled from 
the valence band. The edge of the conduction band re- 



sponds only to the hydrostatic strain, the total confine- 
ment potential for the electron now becomes!^ 



V(r e ) = V e , off + U c 

= Ve.off + a C £hyd, 



(11) 



where a c = -7.17 eW2H26]f or QaAs and -2.09x-5.08 eV 
for Ini-^Ga^As is the hydrostatic deformation potential 
for the conduction band, and £hyd = £xx + £yy + £22 
denotes the hydrostatic strain Shyd- 

The total confinement potential for the heavy hole be- 



comes 



.21 



V(r h ) = V h , off + U v 
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2 \^XX I" €yy ^ZZJ 1 



/36 



{^-xx ^-yy) ^dc X yi 



d (j^-xz ^yz) • 



(12) 



here a v = 1 + 0.16x eV, and b = -1.8 - Q.2x eV for 
In 1 _ 2 .Ga a; As are the deformation potentials of the valence 
band and d = -3.6 - 1.2a; eV^m 

The result of U c and U v in the half section of Fig. [T] 
are shown in Fig. [IT] We take the opposite value of 
the strain induced hole confinement potential (—U v ) here 
(the heavy hole can be treated as confined in a po- 
tential well as electron). Both figures of U c and — U v 
show that the compressive strain (which is the case for 
Ga^Ini-^As/GaAs quantum ring) makes the potential 
energy of the hole and the electron higher inside the 
ring, so the confinement potentials are smaller than the 
ones without strain. Fig. [TT|a) shows that the strain 
induced potential for the electron is smaller in the area 
with smaller p and z, while for the hole the top area of 
the ring has the lowest strain induced potential. As a 
result, the strain pushes the electron and the hole apart 
from each other. In order to see more clearly the effect 
of the strain on the confinement potential and its distri- 
bution, we plot U c and U v along the direction 1 and 2 as 
shown in Fig. [lj for different values of p and z in Fig. 12 



Figure. [12] shows that, as the In concentration is dif- 
ferent in the z direction, the band offsets (confinement 
potential without strain) for both the electron and the 
hole reach their minimum in the top area of the ring 
where z is larger. But in the plane with the same value 
of z, the confinement potential is the same everywhere 
inside the ring. However, in the presence of strain, the 
spatial variation of the total confinement potential for 
the electron and the hole is quite different. In the radial 
plane the confinement potentials are no longer constant. 
The electron, as shown in the first column of Fig. [T2| 
always has a lower potential energy when p e is smaller. 




P h (nm) 



FIG. 11: (Color online) (a) The strain induced shift of the 
band offset of the electron, U c - (b) The strain induced shift 
of the band offset of the hole, — U v . The data taken are in the 
p-z plane. 



Contrary, as shown in the third column of Fig. [12] the 
minimum potentials of the hole are always in the place 
where ph is around 15 nm. In addition, we can find in 
the second column of Fig. [12] that the total potential of 
the electron increases with increasing z e , whereas, the 
confinement potential for the hole still shows a similar 
dependence on z^ as in the case without strain. As a 
result of the strain, the distribution of the confinement 
potential for the electron and the hole becomes very dif- 
ferent and consequently the single particle wave function 
for the electron and the hole will also be different, which 
will affect the polarization of the exciton. 

With the obtained strain induced potential, we solved 
the single particle Hamiltonian in the three dimensional 
ring and we show the results for the energy in Fig. |13| We 
can clearly see that the electron and hole spectrum shows 
similar patterns as in the case without strain (Figs. 13 (a) 
and (b)), but they are quite different when the strain is 
taken into account. The first two transitions in Figs. [13] 
(a) and (b) do not occur at the same B-value, but the 
difference is very small, especially by comparing them 
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FIG. 12: The first (second) column are the electron confinement potential in the p (z) direction for different values of z (p), and 
the last two columns are the results for the hole. Blue-solid curves are the band offset (confinement potential without strain) 
of the electron or the hole, while the red-dashed curves are those for the case when strain is included. Here the confinement 
potentials are in units of eV. 
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FIG. 13: (Color online) Sing particle ground state energy 
for different values of angular momentum as a function of 
magnetic field, (a) and (b) are the six lowest electron and hole 
energy levels without strain, while in (c) and (d) strain was 
taken into account. The arrows indicate angular momentum 
transitions in the ground state. 
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FIG. 14: (Color online) Exciton ground state energy as a 
function of the magnetic field without strain (blue dash line, 
with j/-axis labeling on the left) and with strain (green solid 
line, with y-axis labeling on the right). Inset: the second 
derivative of the exciton ground state energy with respect to 
the magnetic field. 



to the bottom two figures. As shown in Figs. 13 (c) 
and (d) the first transition for the electron takes place 
for B around 7 T, while for the hole it takes place for 
B around 3 T. Moreover, we only have one transition 
for the electron within B = 15 T, the electron is in the 
region with very small radius while the hole prefers the 



top of the ring, as we forecasted. Thus, as a result of 
the large strain distribution difference, the wave function 
distributions for the electron and the hole are much more 
different when strain is present. 

Since the ring is small that has large confinement, the 
Coulomb interaction energy is smaller than the kinetic 
energy, and the exciton will be more polarized as is veri- 
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fied in Fig. 14 where the total exciton ground state en- 
ergy with and without strain are shown. The ground 
state energy is parabolic like with increasing magnetic 
field in the absence of strain, and the amplitude of the 
AB oscillation is small. When strain is included, the 
exciton ground state energy is no longer parabolic-likc. 
The AB oscillation is more pronounced as seen from the 
inset of Fig. [Mj where the second derivative of the exci- 
ton ground state energies with respect to magnetic field 
is shown. By taking the strain into account, the AB 
oscillation is obviously enhanced, but the period of the 
oscillation becomes larger, as the strain induced poten- 
tial confines the electron more towards the center of the 
ring. 



In Fig 15 we show the electron and hole effective posi- 
tion difference for both cases with and without strain. It 
is clear that < Z] x > — < z e >, in case strain is included, 
is more than 5 times larger than the one without strain, 
moreover, < r/j > — < r e > is almost 6 times larger in the 
case strain is included. Thus strain polarizes the exciton. 
Comparing with Fig.[8]of previous section, we notice that 
the electron and hole effective position differences become 
much larger in an InGaAS / GaAs quantum ring, as a re- 
sult of the In concentration difference in the z direction 
and the strain. The electron and hole effective position 
differences in a strained ring even outperforms the re- 
sult of Fig. [8^c), where a strong perpendicular electric 
field is applied. However, similar to the previous section, 
both the electron and hole effective position differences 
decrease with increasing magnetic field, which indicates 
that the exciton decrease its polarity and weakens the 
AB effect. 

When applying a perpendicular electric field, we have a 
similar tunable AB effect as in previous section. The pe- 
riod of the AB oscillation as a function of applied electric 



field is shown in Fig. 16 the dependence of the AB am- 
plitude (we define it by the difference of the value of the 
second derivative of the total energy at B — from the 
value at the magnetic field where the first transition takes 
place) on the electric field is also plotted, but with y-axis 
labeling on the right. The period of the AB oscillation, as 
we found in the case of a GaAs/AlGaAs quantum ring, 
decreases when we increase the top to bottom directed 
electric field. The electron, which determines the period 
of the exciton AB oscillation (see Appendix A), prefers 
to stay close to the center of the ring due to the strong 
strain induced confinement potential. When we apply 
a bottom to top directed electric field, the electron is 
pushed much closer to the center with a smaller effective 
radius, thus the period increases; and as the electric field 
pushes the hole in the opposite direction, the enlarged 
polarity of the exciton should increase the amplitude of 
the AB oscillation, which can be seen from Fig. [THJ If 
we change the direction of the electric field, the period 
will decrease and the AB amplitude will weaken, since 
the electron and the hole are pushed towards each other. 
However, the AB amplitude increases when the electric 
field is so strong (here, when larger than 150 kV/cm) that 
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FIG. 15: (Color online) (a) < Vh > — < r e > (solid line, y- 
axis labeling on the left) and < Zh > — < z e > (dashed line, 
t/-axis labeling on the right) as a function of magnetic field B 
in case strain is excluded, (b) The same as (a) but including 
strain. 



the electron and the hole switch their position: the elec- 
tron attains a larger effective radius than the hole. The 
bottom to top directed electric field has a larger effect on 
the AB oscillation, because the strain induced potential 
counteracts (enhance) the influence of the top (bottom) 
to bottom (top) directed electric field. When comparing 
with Fig. [9j the effect of the electric field here is much 
smaller than in previous section as the height of the ring 
is smaller. 

Our results for Ga x In 1 _ 2 ,As/GaAs quantum rings 
agrees well with recent experimental data on similar 
Ga^Ini-^As/GaAs ringiP^, except for the period of the 
AB effect. The period difference here is because the size 
of the ring in the experiment is a little different from the 
size we took. Figure [T6| shows similar results as Fig. 5 
in Ref. [5S]. by decreasing the top to bottom directed 
electric field, the period of the AB oscillation increases. 
We notice that the second period of the AB oscillation in- 
creases more slowly, and does not take place for magnetic 
field value three times larger than the first one, which is 
very different from the case of an ideal one dimensional 
ringSI. (The explanation for this diffeence is given in Ap- 
pendix A). Furthermore, by studying Fig. (8) of Ref. 
|28j . we know that the PL intensity first increases when 
we increase the top to bottom directed electric field, after 
reaching its maximum value, the PL intensity decrease 
with increasing electric field. As we showed previously, 
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FIG. 16: (Color online) Solid line marked with circle symbols 
is the period of the AB oscillation as a function of applied 
electric field, while solid line marked with prism corresponds 
to the AB amplitude (y-axis labeling on the right). 



has a much smaller effective mass and determines the 
exciton effective radius. As we specified above, this tun- 
able AB effect can only be realized in certain quantum 
rings. We also calculated the exciton energy for a large 
quantum ring in which the Coulomb interaction is larger 
as compared to the confinement energy. In this case, no 
tunable AB effect was found. 

In addition, the strain inside the self-assembled quan- 
tum ring changes the confinement potential of the elec- 
tron and the hole, and makes the confinement potential 
of the electron and the hole quite distinct inside the ring. 
As a result of strain, the polarity of the exciton is in- 
creased and the Aharonov-Bohm effect is enhanced. 

Our results on the tunable AB oscillation of a neutral 
exciton in a quantum ring can be verified experimentally. 
With the help of a perpendicular electric field, it should 
be possible to observe the optical AB effect in small semi- 
conductor quantum rings with relative large height. 



the electron stays in the inner part of the ring while the 
hole prefers the outer part of the ring when the electric 
field is absent, and the exciton has a large polarity. By 
increasing the electric field, the electron and hole move 
towards each other which decreases the polarity of the 
exciton. When the electric field is large enough to switch 
the position of the electron and the hole, the polarity 
of the exciton will increase with increasing electric field. 



This is collaborated by the results shown in Fig. 16 where 
the AB amplitude increases at high electric field the PL 
intensity in Ref. [25] decreases. 



V. CONCLUSIONS 

In this paper we calculated the single particle energy 
of a semiconductor quantum ring and found that the AB 
oscillations of the single particle energies can be tuned 
by an applied perpendicular electric field when the ring 
dimensions satisfy some constraint conditions.. 

For the neutral exciton, we did not find any oscillation 
in the ground state energy of the exciton for both small 
and large quantum rings, as found previously in Refs. 4 
and 5. The AB oscillation can only be seen in the second 
derivative of the ground state energy. But in the presence 
of a strong perpendicular electric field, for small quantum 
ring sizes with a large height in which the effect of the 
electric field is large, the ground state energy of the exci- 
ton shows a weak AB effect as the electric field polarizes 
the neutral exciton in the vertical direction. Moreover, 
this oscillation can be enhanced and the oscillation pe- 
riod can also be changed by increasing the electric field, 
which results in a tunable AB effect. The top to bot- 
tom directed electric field decreases the period of the AB 
oscillation, but has a smaller influence on the oscillation 
than the bottom to top directed electron field. We also 
found that the oscillation of the exciton ground state en- 
ergy mainly originates from the electron motion, which 
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Appendix A: Dependence of the period of the AB 
oscillation on the electron and hole effective radius 

In this appendix, we prove that for the ground state 
with total angular momentum L — 0, the period of the 
AB oscillation is more related to the effective radius of 
the electron, and the effect of the electric field on the 
exciton AB effect is almost the same as on the electron. 

The total Hamiltonian of the exciton is: 
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} + V t . (Al) 



As the In concentration is different in the z direction, m t 
is a function of z. V t in Eq. (Al) stands for the total 



confinement potential which includes band offset, strain 
and so on, and U is the single particle angular momentum. 

In our configuration integral calculation of the exciton 
energy, the l e — and lh — basis function contributes 
most to the L — state before the first momentum tran- 
sition, and the total energy of the exciton is: 



Eq — E, 



B 2 e 2 p\ 



e,0 



kin 
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B e Ph,0 

8m; 



+ v t , 



(A2) 



where E^ in is the total kinetic energy, and p e ^ (ph,o) is 
the effective radius of the l e — (lh — 0) state. Here the 
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effective radius is used, the totai energy is approximate, 
but the difference with the exact resuit turns out to be 
small and we are able to estimate the effect of the electric 
field correctly. 

In the presence of an electric field F — F a , assume the 
effective radius of l e = 0, —1, —2 [lh = 0, 1, 2) to be p e ,o a , 
Pe,u and Pe,2a (f>h.0a, Ph,ia and p h ,2a)- Since the terms 
with B 2 in Eq. (A2| are much smaller (~ 10~ 3 — 10~ 4 ) 



than the other terms, we neglect them here (however, 
they are fully considered in the numerical calculation of 
the total exciton energy). 

At F = F a , the first transition will happen for 
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heB heB 
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2m P 



2m h ' 



2m eP e ,la 2m hP h A a 

(A3) 



and let C Xa = E kmfl - E kinA + V t ,o - V t ,i, as the state 
with larger angular momentum has larger p, the result 
of the first two terms in C\ a is larger than zero, but as 
Vtfl — Vt,\ could be smaller than zero, the sign of C\ a 
is uncertain, however, this term does not have a large 
effect on the transition point. From Eq. (A3) we obtain 



the value of the magnetic field where the first transition 
takes place: 
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When the electric field is F a , the second transition 
takes place when 
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which results in 
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where we neglect the kinetic energies and the potential 
terms. From Eq. (A8) we see that the second transition 
is also dominated by the electron effective radius, and as 
P2a is larger than pi a , the second transition point could 
be smaller by a factor 3 as compared to the first one 
(while in a one dimensional ring, p 2a = P\a)i t ne second 
transition takes place at a magnetic value which is three 
times larger than the first one. 

The above analysis is only a simple estimate, the real 
case will be much more complicated, as a result of the 
coupling of states with different angular momentum. The 
found estimate for the magnetic field transition is approx- 
imate but it turns out to be close to the exact one when 
the particle is more confined. 



here A m — e/H. As the effective mass of the hole rrih is 
much larger (around 9 times) than the electron m e and 
the effective radius of the hole is larger than the electron, 
we see from Eq. ( |A4[ ) that the first transition point is 
determined by the effective radius of the electron. 

Similarly, the first transition point at F = Fb (Fb > 
F a ) should happen for 
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Then the difference of the first transition point by in- 
creasing the electric field is: 



AB 1 = B la - B lb w 1 r x 

A m (m e + m h ) 



Pe.lb Pe,la 
, PeA.aPe.Xb 



Ph.,lb Ph,la 
PhAaPhAb 



(A6) 



Taking into account the effective mass and radius of the 
electron and the hole, it is easy to deduce that the change 
of the first transition is determined by the change of the 
electron effective radius, although the change of the hole 
effective radius may be comparable or even larger than 
the one of the electron (as shown in the previous section). 
That is why the magnetic field of the first transition point 
decreases with increasing electric field, although the ef- 
fective radius of the hole decrease. 



Appendix B: Strain calculation using Eshelby's 
theory 

Eshelby introduced the theory of inclusion in his pa- 
per of 1957^, where he calculated the stress distribution 
for an isotropic and homogeneous inclusion in an infi- 
nite isotropic body. The total strain acting upon the 
inclusion consists of the mutual strains exerted by the 
inclusion and the surrounding material. Here a volcano 
like Ini-^Ga^As quantum ring is embedded in an infinite 
lattice of GaAs and we proceeded as follows: 1) In an 
infinite sample of GaAs, we extract a volcano like ring 
which has the size of the Ini-^Ga^As ring, and leave 
a cavity in the sample. 2) We transmute the material 
of the volcano like ring from GaAs to In 1 _ a; Ga :I .As in a 
strain-free environment. As the lattice constant of GaAs 
is smaller, the ring expands by a fraction cq, which is 
given by the lattice mismatch between the two material, 
i.e. e = ain^xGaxAs/acaAs - 1. 3) By applying hy- 
drostatic pressure, the volume of the ring is reduced to 
its original value. 4) The compressed Ini-^Ga^As is put 
back into the cavity and we allow it to relax. This causes 
radial displacements in both the ring and the surrounding 
material. 

From the point of view of Eshelby's theory of inclu- 
sion, our volcano like ring is the inclusion, embedded in 
a stress-free barrier and characterized by an eigenstrain 
which is material dependent. When the ring is removed 
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from the barrier material, it expands due to this eigen- 
strain. In this case, this eigenstrain is nothing else than 
the lattice mismatch. For small deformations, the to- 
tal strain efj is the sum of the elastic strain and the 
non-elastic eigenstrain e*^: 

(Bl) 

and the total strain is related to the displacement 
through the linearized kinematic relations^ 



1 



(v,i 



i)), 



(B2) 



where Ui is the field of displacements, to be calculated. 

The elastic strain is related to the stress au through 
Hooke's law, so 

a ij = Qjkl e kl = Cijkli^kl ~ e *kl) 

= -^Cijki(uk,i + ui t k) — Cijki^kh (B3) 

in which Cijki denotes the stiffness matrix which is sup- 
posed to be constant!^. When the system is in equilib- 
rium and no external forces are present, we have 

9(Xi i - n *± 1 r ( duk ' 1 m dui ' k \ - r de * kl cms 



dx,; 



dxi 



As a simplication, we restrict ourselves to the case when 
only the change of volume is important, thus we have 
e* kl = e (5fei inside the ring(inclusion) and zero in t he b ar- 



ner. As a result, the right hand side term in Eq. (B4) is 
always zero except at the boundary of the ring. 

Furthermore, we integrate Eq. (B4| from a point r% 



close to the boundary f2 in the ring to a point T2 in the 
barrier (rij below is the normal in the related direction), 



T2 1 r du kt du Lk . 

2 ijkl (ihr ^^> dx i nj = 



dxj 



r)f* 

C , , kl dx n 
dx ■ J 



which after integration becomes: 



(B5) 



\ n i y(Cijki{uk,i + ui,k)) n - (C ljk i(u k .i + ui, k )) r ^j 

= C ijkl e S kl nj(B6) 



Then by solving the 3D equation 
• (CVu) = 



(B7) 



numerically together with the boundary condition ( B6 1 , 
we obtain the displacement u, and the strain efj. 
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